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Abstract 



Singularities in the dark energy late universe are dis- 
cussed, under the assumption that the Lagrangian con- 
tains the Einstein term R plus a modified gravity term of 
the form R a , where a is a constant. It is found, similarly 
as in the case of pure Einstein gravity [I. Brevik and O. 
Gorbunova, Gen. Rel. Grav. 37, 2039 (2005)], that the 
fluid can pass from the quintessence region (w > —1) 
into the phantom region (to < -1) as a consequence of 
a bulk viscosity varying with time. It becomes necessary 
now, however, to allow for a two-fluid model, since the 
viscosities for the two components vary differently with 
time. No scalar fields are needed for the description of 
the passage through the phantom barrier. 
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1 Introduction 

The possibility of crossing the w = —1 barrier in dark energy 
physics is a topic that has attracted a great deal of interest. It is 
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usually assumed that the equation of state for the cosmic fluid 
can be written in the form 



Here the equation of state parameter w = — 1 or p = —p corre- 
sponds to a "vacuum fluid", with bizarre thermodynamic prop- 
erties such as possibly negative entropies PQ. Cosmological ob- 
servations such as SNe la 2J, WMAP 0, SDSS g] and the 
Chandra X-ray Observatory indicate that the present uni- 
verse is accelerating. Moreover, based upon the observed data 
it has been conjectured that w is a varying function of time. 
For instance, as discussed in Ref. jSj, w might have been around 
at redshift z ~ 1 and may be around -1.2 today. (Cf. also 
the very recent analysis of observational constraints on the dark 
energy [7j.) Perhaps is w even an oscillating function in time. 
In view of these circumstances the analysis of a possible crossing 
of the phantom barrier w — — 1, from the quintessence region 
(— 1 < w < —1/3) into the phantom region w < — 1, is obviously 
of physical interest. It ought to be noted that both quintessence 
and phantom fluids lead to the inequality p+3p < 0, thus break- 
ing the strong energy condition. 

The occurrence of a phantom fluid leads inevitably to a sin- 
gularity in the future, called the Big Rip [SJ El HD1 HI] • Often, 
the crossing of the phantom barrier is described in terms of 
one or more scalar fields, besides the gravitational field. For 
instance, as shown recently in Ref. [T2], one can consider in- 
flation, dark energy and dark matter under the same standard: 
phantom-nonphantom transitions may appear such that the uni- 
verse could have been effectively quantum-like both at early 
times and will behave similarly at late times. Cf. also the re- 
lated papers [7J HH]. Here, only one scalar field is considered. 
In several papers - cf., for instance, Refs. [HI E3 IH3 UH HE] - 
cosmological models with two scalar fields (or fluids) are con- 
sidered. In Ref. jH], arguments based upon a stability analysis 
are given why phantom barrier crossing cannot be effectuated 
in terms of one scalar field alone. 

One interesting alternative to a scalar field theory is to allow 
for some sort of modification of the standard Einstein theory. 



p = wp = ( 7 - l) p . 
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One kind of approach is to introduce higher order derivatives of 
scalar fields in the Lagrangian |2D] • The form of modified 
gravity that we shall be concerned with here is to introduce a 
term R a in the gravitational action, where a is a constant. A 
recent introduction to modified gravity of this kind (including 
also Gauss-Bonnet gravity) is given by Nojiri and Odintsov |21| . 
Some other papers dealing with modified gravity from different 
viewpoints are Refs. [221 123 121123 123 123 • 

The body of literature in this field of research is large, and 
the above list of references is not intended to be exhaustive. We 
turn now, however, to the main theme of the present paper, 
which is to introduce a bulk viscosity ( in the cosmic fluid. Bulk 
viscosity is compatible with a universe that is spatially isotropic. 
We shall take ( to be dependent on time. In the paper [2H] 
dealing with Einstein's gravity, we showed how the ansatz of 
letting ( be proportional to the scalar expansion 9 = 3H can 
drive the fluid into the phantom region even if it starts from 
the quintessence region in the non-viscous case. In Ref. j2H] we 
considered a more general situation with the action 

S= ^j d ' X ^ (foR a + An), (2) 

where k 2 = 8irG, /o and a are constants, and L m is the matter 
Lagrangian. This is the kind of action recently analyzed by Ab- 
dalla et al. j3U] • Based upon a natural assumption for the time 
variation of the scale factor, we verified that the field equations 
are satisfied for general a, and we investigated the ansatz of let- 
ting C, be proportional to Q 2a ~ x . This kind of modified gravity 
was further considered in Ref. |31j . We found the same kind 
of behavior also in this case: the fluid can in principle make 
a viscosity-generated passage from the quintessence region into 
the phantom region. There are some obvious advantages of this 
kind of theory: the theory is mathematically simple, it makes 
use of standard concepts from fluid mechanics only (in addition 
to gravity), and there is no need of introducing a scalar field. 

The approach of Ref. [31] is however physically incomplete, 
in the following sense: it involves the gravitational correction 
term R a only, in the Lagrangian. It would be more appropriate 
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to include the Einstein term R also, so that the action takes the 
form 

S = J^j d4x V=g(R+foR a + L m ). (3) 

This form of action is the starting point for the investigation in 
the present paper. As we shall see, the fluid possesses essentially 
the same property as before; it can slide through the phantom 
barrier into the phantom region as a consequence of viscosity. 
But there is one essential difference: the theory now requires 
a two- component model of the fluid. The two components be- 
longing to the Einstein gravity and to the modified gravity vary 
differently with time. It becomes natural to speculate if this 
necessity of dividing the cosmic fluid into two components is 
not after all another demonstration of the difficulties with a 
one-component model that have been found previously in other 
investigations, within a conventional scalar field approach. 

The present series of works is of course not the first time that 
the viscosity concept has been introduced in cosmology. Misner 
was probably the first to introduce viscosity in the cosmic fluid 
in his study of how initial anisotropics in the early universe be- 
come relaxed [32] • Another early paper is that of Padmanabhan 
and Chitre [33]. An extensive review of the development up to 
1990 is given by Gr0n [33]. We have ourselves dealt earlier with 
viscous entropy production in the early universe [33], and with 
viscous fluids on the Randall- Sundrum branes j3J3 EH]. Ren 
and Meng recently discussed a cosmological model with a dark 
viscous fluid described by an effective equation of state, and 
compared with SNe data [3B]. Cataldo et al. considered vis- 
cous dark energy and phantom evolution within the framework 
of the Eckart theory of relativistic irreversible thermodynamics 
|39j . Kofinas et al. considered the crossing of the w = —1 bar- 
rier from a brane-bulk energy exchange scenario containing an 
induced gravity curvature correction term |3D| . As discussed by 
Nojiri and Odintsov ^3] and by Capozziello et al. [7j, a dark 
fluid with a time dependent bulk viscosity can be considered as 
a fluid with an inhomogeneous equation of state. 

In the next section we present the field equations, and the 
conservation equation for energy, and derive herefrom two ex- 
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pressions, Eqs. (fTtjj) and (JT7J), for the coefficient B defined in 
Eq. (f"H?|). A positive value of B leads to a Big Rip. Consis- 
tency of the formalism requires that the coefficients te and r a 
introduced in Eqs. (fT^j) and ()15|) are related to each other. Ex- 
amples with positive values of a are discussed in Sect. 3. The 
case of negative a, discussed in Sect. 4, is exceptional, since the 
physical quantities associated with the modified fluid compo- 
nent does not go to infinity at Big Rip, but rather to zero. Still, 
the Hubble factor, as well as the physical quantities associated 
with the Einstein component, diverge. 

2 General formalism 

We assume the spatially flat FRW metric 

ds 2 = -dt 2 + a 2 (t)dx 2 , (4) 

put the cosmological constant A equal to zero, and adopt the 
energy-momentum tensor of the viscous fluid in the standard 
form 

Tf XV = pU^Uv+ph^, (5) 

where = g^ u + U^Uy is the projection tensor and p = p — 
(6 the effective pressure. In comoving coordinates, the four- 
velocity of the fluid is U° — 1, U l — 0. From variation of the 
action we obtain the equations of motion 

~g^R{l + foR"' 1 ) + RA 1 + afoR 01 - 1 ) 

-afoV.VvR"- 1 + af g^V 2 R a ^ = k% u , (6) 

where T^ u corresponds to the term L m in the Lagrangian. 

Of main interest is the (00)- component of this equation. 
Using that 

R Q0 = -3a/a, R = 6(H + 2H 2 ), (7) 
as well as T 00 = p, we obtain 

3H 2 + l -f R a - 3af (H + H^R^ 1 
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+3a{a - l)foHR a ~ 2 R = K 2 p. (8) 

An important property of Eq. (JBJ) is that the covariant diver- 
gence of the LHS is equal to zero [JT]. Thus 

V%„ = 0, (9) 

just as in the case of Einstein's gravity. Energy-momentum 
conservation is a consequence of the field equations. Contracting 
the expression (j5J) for the energy-momentum tensor with we 
obtain the energy conservation equation 

p+(p + p)3H = %H 2 . (10) 

We now differentiate the expression (jSJ) with respect to t, and 
insert p from Eq. (fTUj) . After some calculation we then obtain 

6H + 9 7 # 2 + ^lf R a - 3a/ [(3 7 - 2)H + 3 1 H 2 ]R a ~ 1 



+3a(a-l)f [(3^-l)HR+R}R a - 2 +3a(a-l)(a-2)f R 2 R a - 3 = Qk 2 C,H. 

(11) 

Recalling that R = 6(H + 2H 2 ), we see that this is a nonlinear 
differential equation for H(t). We shall seek for solutions that 
are compatible with the following basic ansatz for the Hubble 
parameter: 

TT 

H=^-, where X = l-BH t. (12) 

Here B is a nondimensional quantity whose value is dependent 
on a. We take the initial time to be to — 0, and give a subscript 
zero to quantities referring to this instant. If Big Rip shall occur 
(H — > oo), B has to be positive. 

From the mathematical structure of Eq. (jllj) it is clear that 
we cannot eliminate the time dependent terms in the governing 
equation for B simply by letting Q be proportional to a power 
of the scalar expansion 9. This feature contrasts that found in 
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earlier papers j2H] and (31]. It becomes now necessary to allow 
for a two- fluid model. We shall write £ = ((t) as a sum of a 
term referring to Einstein gravity and a term ( a referring to 
modified gravity: 

( = (E + (ai (13) 

where 

(e = r E 6 = 3t e H, (14) 

(a = rj 2 *- 1 = r a (3H) 2a -\ (15) 

te and r a being constants. Then Eq. can be satisfied for 
the Einstein part and the modified part separately. We get 
in this way two different algebraic expressions determining the 
constant B: 

B = -^ + ^k 2 t e , (16) 
(B + 2) a " 1 |9(2 - 0)7 + 3[a + 37 + a(2a - 3) (37 - 

, „i 18k 2 /3\ a , N 

+6a(a-l)(2a-l)i? 2 ) = -^f-J r Q ; (17) 

the time-dependent terms drop out. If /o = (Einstein gravity 
only), Eq. (|16|) is in accordance with Ref. [28J. If the Einstein 
term is absent, Eq. (fTTj) agrees with Ref. pTTj . 

One important property follows at once from the compatibil- 
ity of Eqs. (|T5j) and (JT7J): there must be a relationship between 
t e and T a . This relationship r a = r a (r E ) depends on the values 
of a and f , as well as on the fluid parameter 7 = w + 1. There 
seems to be no direct physical reason behind this dependence; it 
results simply from the mathematical consistency of the model. 

In the next section we shall turn to considering specific mod- 
els. We note, however, the general expressions for B that follow 
directly from the energy conservation equation (jlOj) : as this 
equation has to be fulfilled for the two fluid components sepa- 
rately, 

Pe + (Pe+Pe)3H = % e H 2 } (18) 
Pa + (p a +p a )3H = 9( a H 2 , (19) 
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we get 



B 



3 7 21r E Hl 
2 2 p 0E ' 



(20) 




(21) 




POE 

X 2 ' 



POa 

X 2al 



(22) 



In order to calculate B from the expressions (J2U|) or (|21jl. at 
least one of the initial densities poe or p 0a have to be known. 

Equation (|2T)|) agrees with Eq. ([THjl in view of the first Fried- 
mann equation 3ifg = k 2 Poe at £ = 0. 

3 Positive values of a: Some exam- 
ples 

3.1 Small deviations from Einstein's gravity 

It is natural, for exemplification, to start with a modified gravity 
system where the exponent a in the Lagrangian (jSJ) is small. We 
shall take a to be positive, and also assume that fo is positive 
and small: 



As approximatively R a = 1 + a In R, we see that terms contain- 
ing f a are of second order and thus negligible. The product 
foR a can thus simply be replaced with f R, and the action (JHJ) 
reduces to 



From this equation, or directly from Eq. it follows that the 
Ricci tensor can be taken to be the same as in Einstein's gravity. 
The energy-momentum tensor becomes however rescaled with 



a < 1, f < 1. 



(23) 




(24) 
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the factor 1/(1 + / ) (or 1 — / ). The field equations © reduce 
to 2 

1 1 + J o 

The first Friedmann equation (corresponding to p = v = 0) now 
becomes 

3# 2 = -jLp. (26) 
1 + J o 

As for the determination of the constant B we obtain for the 
Einstein component the same equation (J16)) as before. For the 
a-component (fT7|). when combined with (fTBj) . we obtain the 
simple relation 

Ta = foT E , (27) 

valid for all values of 7. Note the dimensions: [fo ] = cm 2 ^ 1 ', [t e ] 
cm -2 , [r a ] = cm 2 '"" 2 '. (For clarity we write r a instead of T\.) 
Under the present conditions, r a -C te- 

Already from Eq. ()16j) it is clear that B can be greater than 
zero, thus leading to a Big Rip, if r E is big enough. Let us 
summarize the time dependencies of the physical quantities: 

& = ^, PE = % (28) 

t = ^M£, fc = £* (29 ) 
with X = 1 — BHot as before. Moreover, p = pe + p a for all t. 



3.2 The case a =1/2 
This case, corresponding to 

S = 2^ / (^ + 70^+ L m ) , (30) 

is mathematically simplifying, and can moreover be of physical 
interest. Equation ()17|) yields the following quadratic equation 
for B: 

B 2 + 6(7 - 4T 2 )B + 97 2 - 48T 2 = 0, (31) 
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with 



T = —r a . (32) 
Jo 

The product of the two roots for B is equal to 97 2 — 48T 2 . If 
we put 7 = 0, which is the most interesting case, we thus see 
that the product is always negative and there is one positive 
and one negative root. The positive root necessarily leads to a 
Big Rip. When 7 = we get from Eqs. (|3~Tj) and (JH2J), when 
taking Eq. (JT5|) into account, 

9foTE (33) 



4^12 + 27k 2 t e 

Still, r a is of order /oTe, what is physically reasonable. 

In this case and pe vary with time as in Eq. (|28jl . whereas 

( a = r a = const, Pa = ^-. (34) 

The bulk viscosity corresponding to the a- fluid component is 
thus a constant. The density p Q decreases more slowly with time 
than does pe oc X~ 2 . 



3.3 The case a = 2 

This case is quadratic modified gravity, in its simplest form. 
From Eq. ()17|) we get the following cubic equation for B: 

B 3 + (2 + jj 7 ) B 2 + \ lB - 9 -T = 0. (35) 

Let us put 7 = 0. If we draw a curve representing the expression 
(B 3 + 2B 2 — 9T/8) versus B, we see that it has a local maximum 
at B = —4/3 and a local negative minimum at 5 = 0. There is 
thus one single positive root of the equation, for all positive T. 
This root leads in turn to a viscosity-generated Big Rip. A more 
detailed discussion of this case is given in Sect. 4.2 in [IJT]. We 
only note here the time dependencies of the a- fluid component: 

j- f3H \ 3 p 0a 

Ca = T a I — 1 , Pa = — . (36) 
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4 Negative values of a : a = — 1 



If a takes negative values, the situation becomes qualitatively 
different in the following way: Positive values of B will always 
lead to a Big Rip singularity for the Hubble parameter, as H = 
Hq/X —>■ oo at a finite value of t. Similarly, both (e oc X^ 1 and 
Pe oc X~ 2 will diverge. However, both ( a oc X~( 2a_1 ) and p Q oc 
X~ 2a go to zero at the Big Rip. The Einstein fluid component 
and the a-fluid component thus behave quite differently. 
Let us consider a = — 1 as the most typical example: 

S = ^J ^xV^U+^+lS). (37) 

From Eq. ()17|) we get the following equation for B: 

g 13-9 7 + 8T 19T-16T 

2 3+T 4 3+T ' v ; 

T being being given by Eq. (}3*2*|) as before. 

Consider first the vacuum fluid case 7 = 0. The product 
BiB 2 of the two roots is positive, whereas the sum (Bi + B 2 ) of 
them is negative. There is no positive root, thus no Big Rip, in 
this case. The same is true for 7 < 0. The case of general 7 is 
analyzed similarly by observing that B1B2 has the same sign as 
(16T — 97), whereas (B1 + B2) has the same sign as (97 — 3 — 8T). 
On physical grounds, we expect that T 1. It follows that for 
7 < 16T/9 there are two negative roots, whereas for 7 > 16T/9 
there is one positive and one negative root. In the latter case 
there is thus a Big Rip possible, in the special sense explained 
above. Still, there is no need of a scalar field. 



5 Summary. Remarks on Big Rip clas- 
sification 

Summary. The main purpose of this paper has been to investi- 
gate whether the passage through the phantom barrier, from the 
quintessence region (w > —1) into the phantom region (w < —1) 
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can be described as a viscosity-generated phenomenon in the 
general case when the action is as in Eq. As shown ear- 
lier [2H] such a mechanism works well in the case of Einstein's 
gravity, if the bulk viscosity £ oc 9 where 6 is the scalar expan- 
sion. Also, in the more general case of pure modified gravity 
(Lagrangian of the form R a ), the same mechanism was found 
to work if ( oc 6 2a ~ l |29| I31|. what is a natural generalization of 
our basic ansatz. 

The case where the action is as in Eq. Q is more physical 
than those considered in [23 EH] since the Einstein component 
and the modified component are now combined into one La- 
grangian. It turns out that the system shows the same kind 
of behavior as previously: there exists in principle a viscosity- 
driven passage through the barrier w — — 1. It becomes however 
necessary to introduce a two-fluid model, since the bulk viscosi- 
ties for the Einstein component and the modified component 
vary differently with time. Also, physical quantities such as the 
density vary differently for the two fluid components. 

Perhaps is this kind of behavior yet another indication of 
the necessity of introducing a two-component fluid model in the 
late universe, as emphasized, for instance, by Vikman within 
the framework of the scalar field picture jB]. 

It ought to be mentioned again that the the present model 
requires the two coefficients te and r a to be related. This re- 
quirement follows from Eqs. (j!6|) and (fT7|) . There appears to 
be no direct physical reason why this should be so; rather, the 
condition is a consequence of the mathematical consistency of 
the formalism. 

On the Big Rip classification. There are actually several 
variants of the Big Rip phenomenon, and it is of interest to trace 
out to what category the singularities encountered in the present 
paper belong. In Refs. [HI HE UH] two different types of Big Rip 
were discussed. In Ref. [TT] a classification of the known types 
of Big Rip was made and two new types were discovered, so that 
there are according to this four types in all. Let t = t s be the 
instant when Big Rip occurs. From Ref. ^T] we reproduce the 
following classification: 
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• Type I: For t — > t s , a — > oo, p — > oo and |p| — >■ oo 

• Type II: For t — > t s , a — > a s , p — > p s and |p| — > oo 

• Type III: For t — >• t s , a — > a s , p — ► oo and |p| — > oo 

• Type IV: For t — > t s , a — > a s , p — > 0, |p| — > oo and higher 
derivatives of if diverge. 

Here £ s , a s and p s are constants with a s ^ 0. Type I is the 
case one usually associates with the Big Rip concept, emerging 
when w < —1. 

In the present case, when H is given as in Eq. (|12jh the scale 
factor varies with time as 

a{t) = T77? . (39) 

(1 - BH t) 1/B 

Thus, if B > 0, a — > oo when t — > t s = 1/(BH ). It means that 
only Type I is an actual option in our case. If the parameter 
a > in the action (J3J), we have a Type I Big Rip both for the 
Einstein fluid component and the a-component. If a < 0, the 
Einstein component still belongs to the Type I category, but 
the a-component belongs to a new type, different from the ones 
listed above, since ( a — > 0, p a — > 0, p a = wp a — > 0. 
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